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Introduction. 

The present work includes some of the author's 
original researches on the integer solutions of equations 
and linear systems: 

Ts The notion of "general integer solution" of a 


linear equation with two unknowns is extended to linear 


equations with n unknowns and then, to linear systems. 


25 The proprieties of the general integer solution 





are determined (both of a linear and system). 

35 Seven original integer algorithms (two for 
linear equations and five for linear systems) are 
presented. The algorithms are strictly demonstrated and 
an example for each of them is given. These algorithms 


can be easily introduced in a computer. 


INTEGER SOLUTIONS OF LINEAR EQUATIONS 


Definitions and properties of the integer solution of 
linear equations (1.6): 


Let l.e: 
n 
(1) È aix; = b, with not all a: null and b from Z. 


Again, let h e N, and the functions fi : Z > Z, 





i= 1,n. 


Definition 1 
0 
Xi = Xi, i = 1,n, is the particular integer solution 
0 n 0 
of equation (1), if all xi e Z and X aixi = b. 
izl 








Definition 2 








Xi = fi(kı, ... Ki), i = 1,n, is the general integer 


solution of equation (1) if: 





n 

(a) Xo aifi (kı, ... Kn) “b Y(kır --. Ka) e Z", 

i=1 
0 0 
(b) Irrespectively of fi(xi, ... Xn) there is a 
0 0 

particular integer solution for (1) (kı, ... kn) € Z^ so 

0 0 0 
that x: = fi (ki, ... kp) for all i = 1,n. 


We will further see that the general integer solution 


can be expressed by linear functions. 











h 
We consider tor 1 < i < n, the functions di 2; 
oral 
Ci3K3 + di with all Ci, di € 2. 
Definition 3 
A = (cili, is referred to as the matrix associated 
to the general solution of equation (1). 
Definition 4 
The integers ki, 66. ks, 1 < s < h, are independent 


if all the corresponding column vectors of matrix A are 


linearly independent. 


Definition 5 





An integer solution is s — times undetermined if the 
maximal number of independent parameters is s. 
Theorem 1. The general integer solution of equation 


(1) is (n-1) - times undetermined. 


Proof 


We suppose that the particular integer solution is of 


the form: 
r — 
(2) Xi = È i:P. + Vi, 1 = 1,n, with all Uie, vieZ, 
e=1 
Pe = are parameters: from Z, while a £r € — d. 
0 0 
Let (Xı, ... Xn) a general integer solution of equation 





(1) (ve are not interested in the case vhen the equation 





does not accept integer solution). The solution 
E . AQ rs Eee 
Xj = anlj + Xj for go lə ol 
n-li 0 
Xn = -( X aşk; Xn) 
jc 
is (n 15) times undetermined (it can be easily checked 


that the order of the associated matrix is n — 1). Hence, 





there are n 1 undetermined solutions. let, in the 





general case, a solution n - 1 times be undetermined: 
n-1 -. 
Xi = È CijKj + di i = l,n, with all Ci, di e 2. 
j=1 
The case when b = 0. 
n n 
Then X aixi = 0. It follows X aixi = 
izl izl 
n n-1 n n-1 n 
= X a (X Cik: + di ) = Xa: 22 ike > adi 
izl j-1 izl j-1 izl 
= 0 
n 
For ki, = 1 and k; = 0, 1.. jo, it follows È aici; ... 0. 
izl 


Let the homogenous linear system of n equations with n 


unknowns be: 


ei Ste 
2. XiCij = 0, j —-1,n- 1; 
izl 
n 
25 xidi = 0 
i=1 
which, obviously accepts the solution xi = ai, i= 1,n, 


different from the trivial one. Hence the determinant of 





the system is zero, i.e., the vectors Cj = (cu, ... Cu)”, 
j= 1,n- 1, D= (dı, ..., Gs) are linearly dependent. 
But the solution being n - 1 times undetermined it follows 


that C;, j = I,n - 1, are linearly independent. Then (Ci, 
Cn-1) determines a free sub-module Z of the order n 1 





in Z, of solutions for the given equation. 


Let us see what can be obtained from (2). We have: 
n n r 

0 = aixi = x ai (X uiePe + vi). As above, 
i=1 izl e=1 


we obtain: 


n 
>: aivi = 0 and È iie = 0; and similarly, the 


vectors Un = (um, ., Unn) are linearly independent, 


ha Tti Eo “Uk, h = 1,r, and V = (vz, ., Vn) are 
particular integer solutions of the homogenous linear 


equation. 


Subcase (al) 





Un, h = 1,r, are linearly dependent. It gives 


(Ui, ..., Ur} = the free submodule of order r in z" of 





solutions of the equation. Hence, there are solutions 


from (Vi, ..., Vn-1} which are not from {U1, eş Ur), Chis 
contradicts the fact that (2) is the general integer 


solution. 


Subcase (a2) 





Un, h = 1,r, V are linearly independent. Then, 
(Ui, ..., Ur} + Vis a linear variety of dimension < n - 1 
= dim (Vi, ..., Vn-1} and the conclusion can be similarly 


drawn. The case when b #0. 


n n n-1 
So, È aix; = b Then X a: ( X cişk; + di) = 
i-1 i-1 el 
n-li n n 
= x (E 81€13) Ks; + > aidi = b, V(Ki, ..., Kna, 
j=1 i=1 i=1 
n 
Hen" As in the previous case we get X aidi = b and 
izl 
7 e tt 
Xo aci = 0, V3 = 1I,n - 1. The vectors C, = 
izl 
(Ci, «++, Cnj) , 3 = 1,n — 1, are linearly independent 
because the solution is n - 1 times undetermined. 
Conversely, if Cı, ... Cna, D (where D = İdi; sə da) *) 
n-1 
were linearly dependent, it would mean that D = È sC,, 


n 
with all s: scalar; it would also mean that b= X aid; 


izl 
n n—1 n-1 n 
x. a(X sı) = È s3( Xoaicı) = 0. Impossibly, 
isf j=1 j=1 i=] 
(3) then (Cı, ..., Cn} + D is a linear variety. 
Let us see what we can obtain from (2). We have: 
n n r r n 
b= È aix; = Xa: ( È UiePe +t Vi) = X ( DWaittie) Pe 
i-1 i-1 e=1 e-] izl 
n 
+ 2; aivi 
izl 
n n ——— 
and, similarly: Yaivi = band È aiuie =0, v e = 1,r, 
izl que 
respectively. The vectors Ue = (nj xo; Une)", e = 





1,r, are linearly independent because the solution is r 
— times undetermined. 
A procedure like that applied in (3) gives that Ui, 
., Ur, V are linearly independent, where V = (vı, ..., 
ub “Then. fU ... Ur} + V = a linear variety = free 
sub-module of order r < n- 1. That is, ve can find 
vectors from Cı, ..., Cn- + D which are not from (U,ş, 


., Ur} + V, contradicting the "general" characteristic 


of the integer number solution. Hence, the general 





integer solution is n - 1 times undetermined. 
Theorem 2. The general integer solution of the 
n 
homogenous linear equation È aix; = 0 (all a e ZN (0)) 
izl 


can be written under the form: 


n-1 ——— 
(4) Xi = > Cişki, i= l,n, (with Or 5 o di E 0) . 
j=1 
Definition 6. This is called the standard form of 





the general integer solution of a homogenous linear 
equation. 


Proof 





We consider the general integer solution under the 


n-1 — 
form: x” X jP; + di, i = 1,n, with not all d; = 
ge 
0. We show that it can be written under the form (4). 
The homogenous equation admits the trivial solution 


xi = 0, for 
ə 0 0 n-1 0 
i = 1,n. There is (Pa; Se Pha) gez so that È Cijp; + 
j=1 


tdi = 0, V i = 1,n. Substituting: P; = k: + p;, 


j = I,n - 1, in the form from the beginning of the 
proof we obtain form (4). We have to 


mention that the substitution does not diminish the 


degree of ii: TEE 
generality as P; e Z because kj e 2, j e 1,n - 1. 

Theorem 3. The general integer solution of a non- 
homogeneous linear equation of its associated homogenous 
linear equation + any particular integer solution of the 


non-homogeneous linear equation. 


Proof 
n-1 ——— 


Let xi = X cişki, i = 1,n, be the general integer 


solution of the associated homogenous linear equation 
and, 


again, let xi = vi, i= 1,n, be a particular integer 
solution 
of the non-homogeneous linear equation. Then, 
n-1 
Xi = È Cijky + 
j=1 
+ Vi, i = 1,n, is the general integer solution of the 
n 


non-homogeneous linear equation. Actually, X aixi = 


n-1 n-1 n 
x ai ( X cişki) + XM aiVi 
=1 el i=1 


I 
IM 5 
o 
b 
"M 
Q 
E 

. 
+ 
Le 
$ 
| 


1 el nb 


if xi = xi, i = I,n, is a particular integer solution of 


10 


the non-homogeneous linear equation, then xi = Xi — Vi, 


i = 1,n, is a particular integer solution of the 
homogenous linear 


0 0 
n-1 





equation; hence, there is(kı, ..., kai) € Z so that 
n-l o 0 aun n-1 9 : 
X cişkj = Xi — Vi, Vi = 1,n, i.e., X cak: + vi = xi,Vi € 
gel dei 
1,n, vhat had to be proven. 
n-1 m 
Theorem 4. If xi = X cak, i = 1,n, is the 
j-1 


general integer solution of a homogenous linear equation 
(Ci, ep Cnj) = 1, 


The demonstration is made by reductio ad absurdum. 
If 1 jo, Los ol, 


so that (C15, 1 Cni, ) = di, + t1, 
then Cis, = Gi, dis, with (Cis, px AUS Cn3, ) mE Vi -1,n. 
But xi = Ci; , i = l,n, represents a particular integer 
0 
n n ' n 
Solution as È ax = X aici; = l/dj ! X ic 0 
è ] 0 0 : 0 
izl izl izl 
(because xi = Ci: , i = l,n, is a particular integer 
0 


solution from the general integer solution by 
introducing k; = 1 and 
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0 


k = 0, j # jo). But the particular integer solution 


Xi = 5 i — 1,n, cannot be obtained, by introducing 


vhole number parameters (as it should), from the general 


integer solution, as, from the linear system of n 


equations and n - 1 unknowns, which is compatible, we 
obtain: 
n ' ' —— 
Xi = 2 CijKj + Cin di, ks = C13, po ok = l,n, 
azı 
1220 


Leaving aside the last equation--which is a linear 








combination of the other n 1 equations, a Kramerian 
system is obtained. It follows up: 
C11 Cij ...Cı1,n-1 | 
Cn-1,1- 8 -— e e e Cn-1,n-1 | 
T 
kj = = € 2 
0 ; dj 
0 


Cn-1,1+Cn-1, a di e Cn-1,n-1 | 


The assumption is false and thus the proof results. 


Theorem 5. Considering the equation (1) with 
(a1, ..., an) = 1, b = 0, and the general integer 
solution 
n-1 ——— 
Xi = 2 Cişki, i= l,n, then (ai, ep Qi-1; Aitly e c er 
qr 
An) — (Ci, ..., Cin-ı), Vi O 1,n. The demonstration is 


made by 
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double divisibility. Let io, 1 € iş < n, be arbitrary 
but 


n-1 

fixed. xi X ci,k,. Consider the equation È aiXi 
j=1 j Mi Lu 

= “ai Xi o. We have shown that x: = cı, i = 1,n is a 


particular integer solution irrespective of j, a < j < 


<n- 1. The equation Ó axı = -ai ci j, obviously, 
dəndə ne 0 

admits the integer solution xı = Ci;, i .. iç. Then 
(a1, ..., ai-cı, aisı,..., ön) divides -ai ci Jj, but, as 
we 

0 0 0 0 
have assumed that (ai, ..., an) - 1, it follows that 
(ai, ..., ayı, ai ır 007 an) | eis; irrespective of j. 
Hence (ai, ..., aia, Aisir seer an) | (cı ı, seer Cimi), 
Vi 0 1,n, and the divisibility in one sense vas 
proven. 


Inverse Divisibility 





Let us suppose the contrary and say that > ii 0 


O 1,n for which (ai, ..., ai-cı, Bi x, ..., an) — Gia w= di 
em 
1 1 1 1 


(eiu, ..., Cina); we have considered dii and di 2 
1 1 1 


without restricting the generality. di ildi 2 according 


13 


to 
T 1 


the first part of the demonstration. Hence, >» d 0 Z so 
that 


di 2 =d ! di i, lal ... 1. 
1 1 
n-1 n-1 ' n 
Xi O Ci 4k; d 1 d: ï O Ci ikay O diXi 0 YO 
ak = 
j=1 2 ^ el : i-1 Jal 
n-1 ' 
= -ai Xi O diXi = —dâi d 1 d: 1 O ca jK5; where 
Hs l s : 1 al E il. 
Tl j=1 
(Ci lr , Ci n-1) D 1 
1 ib 
The nonhomogeneous linear equation Ó aixi = -ai di 
admits Acla E. 


integer solution because ai di is divisible by 


0 


(Ai, ..., Bi ip, dişi, «++, An). Let Xi = Xi, i .. ij be 
its 
1 əl: 

particular integer solution. It follows that the 
equation 

n 0 

O aixi = Ü admits the particular solution xi = xi, 

izl 

J li, Xi = di , which is written as (5). We show that 


ib 


14 


15 


cannot be obtained from the general solution by integer 


number parameters: 


n-1 0 
Ó Cişkj = Xi, 1 Ia? 
j=1 
nk 
d ! dia O Ci3k; = di: (6) 
1 ^E 1 
yel 


But equation (6) does not admit integer solution because 


a 1 dia dia thus, contradicting, thus, the "general" 


characteristic of the integer solution. 


As a conclusion we can write: 


Theorem 6. Let the homogenous linear equation be: 
n 
O ax; = 0 with all ai OZ N (0), and (ai, ..., an) - 
Lz 
i=1 
h es 
Let x; = O Cakiı, i = 1,n be, with all ci; 0 Z, all kj 
j=1 


whole integer parameters and h 0 N, a general integer 
solution of the equation. Then, 


1° the solution is n - 1 times undetermined; 





2° Vj ilem. 1 we have qeu e Cn) diş 


0 
3 Vi = 1,n we have (Cu, ..., Cin-1) — (ai, ..., aia, 


16 


Misi, ..., an). The proof results from Theorems 1, 4 and 


Note 1. The only equation of the form (1) which is 


n — times undetermined is the trivial equation 0 ! xı + 
se rəyə səbə xb > ee OL, 
Note 2. The converse of theorem 6 is not true. 
Counterexample: 
X1 = — kı + ka 
X2 = 5ki + 3kə 
X3 = 7ki - ko, ki, k ÜZ (7) 


is not the general integer solution of the equation 


-13x: + 3x; — 4x4 = 0 (8) 





although the solution (7) verifies the points 1°, 2° and 
3 


of theorem 6. (1, 7, 2) is the particular integer 
solution 


of (8) but cannot be obtained by introducing integer 
number 


parameters in (7) because from 


| 
RA 


—kı + ko 


5ki + 3k> = 


| 
J 
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follows that k = 1/2 0 Z and k = 3/2 0 Z (unique roots). 


Reference: 


[1] Smarandache, Florentin, “Whole number solutions of 


linear equations and systems”, M. Sc. Thesis, 1979, 





University of Craiova (under the supervision of 


Assoc. Prof. Dr. Alexandru Dinca). 
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AN INTEGER NUMBER ALGORITHM 


TO SOLVE LINEAR EQUATIONS 


An algorithm is given which ascertains whether a 
linear equation admits integer number solutions or not; if 


it does, the general integer solution is determined. 


Input 
A linear equation aıxı + ... + anXn = b with a: b Ü Z, 
x; being integer number unknowns, i = 1,n and not all 
ai = 0. 
Output 


Decision on the integer solution of this equation; if 
the equation has solutions in Z, its general solution is 


obtained. 


Method 

Step 1. Calculate d = (ai, ..., an). 

Step 2. If d/b then "the equation has integer 
solution"; go on to Step 3. If d/b then "the equation 


does not admit integer solution"; stop. 


Step 3. Consider h: = 1. If |d| .. 1 divides the 


equation by d; 
Step 4. 


that a: = a. 
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consider ai: = a,/d, i= 1,n, b: = b/d. 


Calculate a = min |a,| and determine an i so 


as...0 





Step 5. If a. 1, go on to Step 7. 
Step 6. If a = 1, then: 

(A) Xi = —-(aiXi + ... + arıXıı + anıXısı + 
To əkə a cce ob uar 

(B) Substitute the value of x; in the values of 
the other determined unknowns. 

(C) Substitute integer number parameters for 
all the variables of the unknown values in 
the right term: ki, ko, ..., Kn-2 and Ka, 
respectively. 

(D) Write down the general solution thus 
determined; stop. 

Step 7. Write down all aş, j .. i and b under the 
form: 
aj = aid; + rj 
a: b 
b dıq + r, Where q; = = , q - Seem 
LED tees — ai _ 
Step 8. Write xi = - qiXi — ... — qııXıı — QiıXısı — 


— ... — QnXn + q — th. Substitute the value of x: in the 


values of the other determined unknovns. 
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Step 9. ai: = 4, and ai: = —a: 
ğ b:= r 
Consider . Xi: = th 
: h = h+ il 
ai-l: = r 
Ai : = Tin 
an? = rn and go back to Step 4. 
Lemma 1. The previous algorithm is finite. 
Proof 
Let the initial linear equation be aixi + ... + anXı = 
= b with not all ai = 0; it is considered that min |a,| 
a.s... 
a1 . 1 (if not, it is renumbered). Following the 


algorithm, once ve pass from this initial equation to a 


new equation: aiti + azx2 + ... + anXn = b", with |ai| < 


€ [as| for i = 2,n, |b'| < |b| and aı = an It 
follows that min İa.) € min ļas|. We continue similarly 
a ..0 agal 


s 


and after a finite number of steps ve get, at Step 4, at 
a: = 1 (as, every time, at this step the actual a is 
strictly smaller than the previous a, according to the 
former note) and in this case the algorithm terminates. 


Lemma 2. Let the linear equation be: (25) aıxı + 





+ a2X2 + ... + anX%n = b with min lasl = ai and the 
as...0 
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equation: (26) — aştı + r2X2 + ... + YnXn = r with ti = 
= — Xi — Q2X2 — ... — QaXn + Q where ri = ai — aigi, i = 
c ai b 
—-2,n, r = b — aig while qi = = por m . Then 
— a = — ai şı 
0 0 0 
Xı = Xı, X2 = Xo, ..., Xn = Xn is particular solution of 
0 0 
equation (25) if and only if tı = tı = — xı — qox» Z ... — 
0 0 
QnXn + Q, Xo, ..., Xn = Xy is a particular solution of 


equation (26). 





Proof 
0 0 0 
Xı = Xı, X2 = Xo, ..., Xn = Xn is a particular solution 
0 0 0 0 
of equation (25) 1 aıxı + a>X> + ... + anXn = b ] aixi + 
0 0 0 
+ (r2 + aıqz) xə viven aiqa) Xn = aid + r ] £>X> + 
0 0 0 0 0 
To... T YTaXı ai (-Xi1 -q2X2 eee — GnX¥n + q) = rl — aiti 4 
0 0 0 0 0 
BP BEKO nemi Ke or J). ba Sa xox ee SS Xa.” KASS a 


particular solution of equation (26). 


Lemma 3. xi = Ciiki + ... + Cin-tkn-1 + di, i = 1,n is 
the general solution of equation (25) if and only if: 


(28) 


tı = — (Cı: + d2Cz + ... + GnCni) ki — ... — (Cina + 


+ Q2C2n-1 Tee Ub QnCnn-1) Kn = (dı F qədə F- əə” Sar Inn) + 


+ q, Xi = Cejıkı 5 Cjn-1Kn-1 + dj, j= 2,n 
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is a general solution for equation (26). 


Proof 
0 0 0 0 0 
tı = tı = — Xi — qox? — ... — OnXn + Q, Xə = Xə, ..., 
0 
., Xy = Xn is a particular solution of the equation (25) 


0 0 0 


] Xi = Xi, X2 = Xo, ..., Xn = Xn iS a particular solution of 
0 0 
equation (26) | > kı = kı 0 Z, ..., Kn = kn 0 Z so that 
0 0 0 m. 
Xi € Cika + ... + Cinukni + Gi = Xi, i = 1,n | > k = 
0 0 0 0 
= kı ? Vig. İZ kiık, 0 Z so that xı € Ciki + ... + Cin-iKnu + 
0 eu 0 
+ d; = Xi, a = Zn and tı = — (C4 + qoco c Qe c QnCni) Kı = 
0 
— ... — (Cini + Q2€2a-31 £ ... + GQnCnn-1) Kn-1 — (di + Qəd? + 
0 0 0 0 
tw + qada) doq = -= XQ — Geko — eau — Qaxa t q — ti. 
Lemma 4. The linear equation (29) aıxı + a>X> + ... + 
+ anxn = b with |ai| = 1 has the general solution: (30) 
xı = — (a>k> + ... + ankn — b)aı, 


Proof 


0 0 0 
Let Xı € Xi, X2 = Xo, ..., Xn = Xn be a particular 


0 0 
solution of the equation (29). > k2 = xo, ..., Kn = kn SO 
0 0 0 0 
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that xı: = — (a>X> + + anXn — b)aı = Xi, Xə = Xə, ..., 
0 
ap Mg = Spk 
Lemma 5. Let the linear equation be aıxı + â&>X> + 
+... ax, = b, with min [a = aı and ai = mud i = 2,n. 
as... 


Then, the general solution of the equation is: 


Xı = — (g>k> eee OF Qnkn = q), 


Xi = ki Ü Z, 


Proof 
Dividing the equation by ai the conditions of Lemma 4 
are met. 


Theorem of Correctness. The preceding algorithm 








correctly calculates the general solution of the linear 


equation aıxı + ... + anXn = b with not all a: = 0. 


Proof 
The algorithm is finite according to Lemma 1. The 
correctness of steps 1, 2, and 3 is obvious. At step 4 


there is always min |a;| as not all a: = 0. The 
asas 


correctness of substep 6(A) results from Lemmas 4 and 5, 
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respectively. This algorithm represents a procedure of 
obtaining the general solution of the initial equation by 
means of the general solutions of the linear equation 
obtained after the algorithm was followed several times 
(according to Lemmas 2 and 3); from Lemma 3 it follows 
that to obtain the general solution of an initial linear 
equation is equivalent to calculate the general solution 
of an equation at step 6(A), equations whose general 
solution is given in algorithm (according to Lemmas 4 and 
5). The theorem of correctness has been fully proven. 


Note. At step 4 of the algorithm we consider 


a: = min lasİ so that the number of iterations be as small 
as... 0 
as possible. The algorithm works if we consider 
a: = İail . max İasl but it takes longer. The algorithm 
s-1,n 


can be introduced in the computer. 


Application 








Calculate the integer solution of the equation: 


6X, — 12x; — 8x3 + 22x, = 14 


Solution 


The former algorhythmus is applied. 


L; (6, — 12, - 8, 22) = 2 
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2114 so that the solution of the equation is in 


Z 

h: = 1 

|2| . 1; dividing the equation by 2 we get: 
3xi = 6X2 — 4x3 + 11x, = 7 

a: = min dele. |-4ļ, |11|} =3, i=1 
& ii 

-6 = 3.(-2) + 0 


11 = 3.3 + 2 


7=3.24+ 1 
Xi = 2X2 + 2x3 — 3x, + 2 - ti 
a2: = 0 ai: = —3 
a3: = 2 b: = 1 
a4: = 2 Xi: = tı 
h: = 2 





We have a new equation: 


—3tı + 0 ! Xo + 2x3 + 2x, = 1, 


a = min (1-31, EF |2|} = 2, and 
də 
aw 1 


8. 
the value 
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N 
I 


2. b 4.4 0 


210-40 


HR 
I 


X3 € 2t: - 0 ! X2 — X4 + 0 — to. Substituting 
of x3 in the value determined for xi we get: 





X1 = 2x2 — 5x, + 3t1 — 2t2 + 2 


ai: = 1 a3: = —2 
a2: = 0 b: = 1 
aa: = 0 X3: = t3 

h: = 3 


We have obtained the equation: 


1 ! te + 0x2 — 2t2 + Ül x, = 1, 


a = 1, and 

i = 1 

(A) tı € -(Ü l x? — 2t2 + Ül x, — 1) ! 1- 
2t2 + 1 


(B) Substituting the value of tı in the values 





of xı and x3 previously determined, we get: 


Xı = 2x2 — 5x, + 4tə + 5 and 


X3 = -x4, + 3t2 + 2 
(C) X2: = kı, X4: = ko, t2 = ks, ki, ko, ks 0 2 
(D) The general solution of the initial 


eguation is: 
xı = 2k, - 5k, + 4k; + 5 
X2 = kı 


= -k, + 3k3 + 2 


x 
w 
| 
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kı, ko, ks are parameters 0 Z 
Reference 
[1] Smarandache, Florentin, Whole number solutions of 
equations and systems of equations--diploma paper, 


University of Craiova, 1979. 


28 


ANOTHER INTEGER NUMBER ALGORITHM TO SOLVE 


LINEAR EQUATIONS (USING CONGRUENCY) 


In the present part a new integer number algorithm 
for linear equations is presented. This is more "rapid" 
than W. Sierpinski's presented in [1] in the sense that it 


reaches the general solution after a smaller number of 





iterations. Its correctness will be strictly 


demonstrated. 


INTEGER NUMBER ALGORITHM TO SOLVE 
LINEAR EQUATIONS 


Let us consider the equation (1); 
(the case ai, b are in Q, 


i = 1,n, is reduced to the case (1) by bringing to the 
same denominator and eliminating the denominators). 

Let d = (aij, ..., an). If dib, then the equation does not 
admit integer solutions, vhile if db, the equation admits 
integer solutions (according to a well-known theorem from 


the theory of numbers). 
If the equation accepts solutions and d .. l, we 


divide the equation by d. Then, we can agree that d = 1 
(we do not make any restriction if we consider the maximal 


co-divisor positive). 
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(a) Also, if all a: = 0, the equation is trivial; it 
admits the general integer solution xı = ki Ü Z, i= 1,n 
when b = 0 (the only case when the general integer 
solution is n - times undetermined) and does not have 
solutions when b .. 0. 


(b) If>i, Aad = my so that a; = + 1, then the 


general integer solution is: 


n 
xi = -a; ( Ó ajk; - b) and xs = k, 02, s 0 
j= 
del 
O dle ə SELİ 
The proof of this assertion vas give in [4]. All these 
cases being trivial, we will leave them aside. The 
following algorithm can be written: 
Input 
A linear equation: (2) 
n mx 
Ö aixi = b, ai, b Zp a; .. Rd T = 1,n 
izl 
with not all a: = Ü and (ai, ; aş) = 1 


Output 


The general solution of the equation. 


Method 


h: = 1, p: = 1 
Calculate min (rl, r / ai (mod aş), 
1<i,j<n 
£20 Ir] < las]? 


and determine r and the paid (i,j) for which 


this minimum can be obtained (when there are 





more possibilities we have to choose one of 


them). 
If |r| œ 1 go on to step 4. 
If |r| = 1, then 
n 
SpE (Hasty. =O asxs + b) 
s=1 
só{i, j} 
ai- r n r-a: 
Xj =xr(at+---- ! O asXs + ——— b) 
a; s-1 aj 
só{i, j} 


(A) Substitute the values thus determined of 
these unknowns in all the relations (p), 


p > 1, 2... (if possible). 
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əl 





(B) From the last relation (p) obtained in 
the 
algorhythmus substitute in all relations: 
(p - 1), (p - 2), ..... (1). 
(C) Every relation, starting in order from 
(p — 1) should be applied the same 
procedure as in (B), then (p - 2), ..., 
., (3), respectively. 
(D) Write the values of the unknowns xı, 
i — 1,n from the initial equation vriting 
the corresponding integer number 
parameters 
from the right term of these unknovns 
vith 
kı, Sw; Kad); STOP. 
4. Write the relation (p): 
ai- r 
Xj tn — MS Xi, 
aj 
5i Consider Xj: = th h: = h + 1 
ay: = r p: = p t 1 


The other coefficients 


and variables 





remaining 


unchanged and go back to step 2. 
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The Correctness of the Algorithm 





Let us consider linear equation (2). Under these 
the following proprieties exist: 


conditions, 
(İrl, r / ar (mod aş), 


Lemma 1. The set M — 


0 € |r| € Jaş;şlı has a minimum. 


Proof 





Obviously M d N” and M is finite because the equation 


n and considering 


a finite number of coefficients: 


has 
the possible combinations of these, by twos, there is 


all 


2 
the maximum AR, (arranged with repetition) = n elements. 
Let us show, by reductio ad absurdum, that M.. , 


Me ] ai = 0 (mod aş) V i, j 0 1,n. Hence a; = 
= 0 (mod ai), Vi, j 0 1,n. Or this is possible only 


1,n, which is equivalent to 


when | a; | = lasl, Vi, 230 


(a1, ..., an) = ai, Vi 0 I,n. But (ai, ..., an) = 1 


according to a restriction from the assumption. 


follows that la: = 1,n, Vi 0 1,n a fact which contradicts 


the other restrictions of the assumption. 


it follows that M has a minimum. 


M.. 0 and finite, 
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Lemma 2. If |r| = min M then |r| < lai], Vi 0 
1<i, j<n 


0 1,n.: 


Proof 


We assume, conversely, that »io, 1 € io € n, so that 


| x | > la: E Then | x | > min (lal) - la: | ge 1l, 1 < jo < 
= 0 1«j«n 0 
€ n. Let ap , 1 € po € n so that | 8» | > la: İ and ap is 
0 
not divided by a;. There is such a coefficient in the 


equation as las, | is the minimum and not all the 
coefficients are equal among themselves (conversely, it 
would mean that (ai; sss, an) > ayr = + 1, which is against 


the hypothesis) and, again, of the coefficients whose 


module is greater than | ass, | not all can be divided by a; 


(conversely, it would similarly result that (ai, ..., 
an) = 

= aj; f 1). We write [ap /â5, 1 = qo 0 Z (the whole 
number part), and r = ap, 7 daş, 02. We have ap, / £o 
(mod a; ) and 0 < lro] < laş | < [ai | < |r|. Thus, we 


have found a ro with |ro| < |r| which contradicts the 


definition of minimum given to Ir]. Contrary to the 


assumption. Thus, |r| < Jal, Vi 0 işa, 


Lemma 3. If |r| = min M= 1, for the pair of indices 
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n 
Xi = r (-ajtn - Ó asks + b) 
sel 
só(i, j} 
ai-r n r-ai 
xı =r (aitn + ---- Ó asks + ——— b) 
aj s=1 aj 
só{i, j} 
Ke = ky 0 Zy 6 DAİ; "nb. fig} 


is the general integer solution of equation (2). 


Proof 
m E > ə. 
Let Xe = Xe, e € 1,n be a particular integer solution 
0 
of equation (2). Then» ks »- x, 02, s 0 (1, ..., n) \ 
0 ai- ro 
\ {i,j} and ti = x? 4 ---- xi 0 Z (because a:-r =;aj), so 
aj 
that: 
0 ai-r no 0 
Xi = r -aj (Xj ba SSS Xi) = O asXs + b —Xi 
a: s=1 
só(i,j) 
0 ai-r 0 ai- r No 
Xj; = r -—a:(X: £ ——— xi) iman O asXs + 


35 


só{i, j} 
r-aı 0 
F —— b = x: 
aj 
0 
and x, = ks = x, s 0 (1, ..., n) \ {i,j}. 
Lemma 4. Let |r| .. 1 and (i,j) be the pair of 


indices for which this minimum can be obtained. Again, 


let the system of linear equation be: 


n 
ajtn + rx; + 6 asXs = b 
sel 
s..{i, 3} 
(3) 
ar 
tn EXC oq o X 
aj 
R zx f 
Then, Xe = Xe, e = 1,n is a particular integer solution 


0 


for (2) if and only if x. = Xe, e 0 (1, ., n} \ {j} and 


0 0 ai-r 
th = th = X + -——- X; is the particular integer solution 
aj 
of (3) 
Proof 
Xe = Xe, e = 1,n is a particular integer solution for 
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n 0 Jn 0 0 ai-Yr o 
1 O aexe- b ] O asXs + a4(xj + ———— Xi) + 
e=1 sci a: 
s..{i, Jj} 


] ajtn + rx; + O âsXs = b and ty = Xj + 
s=1 
s..{i, J} 
ai- r 0 
eS xQ 0 E 
aj 
0 
] x. = x, e Ü (1,2, ..., n) \ (1) 
0 
and th = tn is a particular integer solution for (3). 
Lemma 5. The former algorithm is finite. 
Proof 
When |r| = 1 the algorithm stops at step 3. We will 
discuss the case when |r| .. 1. According to the 
definition of r, |r| 0 N*. We show that the row of r - s 


successively obtained by following the algorithm several 


times is strictly decreasing to 1. Let rı: be the first 
obtained by following the algorithm one time. |ril = 1, 


go on to steps 4 and, then 5. According to lemma 2 İr:l < 
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< lail, Vi = 1,n. Now we shall follow the algorithm a 
second time, but this time for an equation in which rı 
(according to step 5) is substituted for ai. Again, 
according to lemma 2, the new |r| written |r;| will have 
the propriety: |r2| € |ri|. We will get to |r| = 1 as 

|r| > 1 and |r| < 4, and if |r| ..1, following the 
algorithm once again we get |r| < |ri|, a.s.o. Hence, the 
algorhythmus has a finite number of repetitions. 


Theorem of Correctness. The former algorithm 





calculates correctly the general integer solution of the 


linear equation (2). 


Proof 
According to lemma 5 the algorithm is finite. From 
lemma 1 it follows that the set M has a minimum, hence 
step 2 of the algorithm has meaning. When |r| = 1 it was 
shown in lemma 3 that step 3 of the algorithm calculates 


correctly the general integer solution of the respective 





equation (the equation that appears at step 3). In lemma 
4 it is shown that if |r| ..1, by the substitutions steps 


4 and 5 introduce in the initial equation the general 
integer solution remains unchanged. That is, we pass from 


the initial equation to a linear system having the same 
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general solution as the initial equation. The variable h 


is a counter of the newly introduced variables which are 





used to successively decompose the system in systems of 
two linear equations. The variable p is a counter of the 
substitutions of variables (the relations, at a given 
moment, between certain variables). 

When the initial equation was decomposed to |r| = 1, 
we have to follow the reverse way: i.e., to compose its 
general integer solution. This reverse way is directed by 
the substeps 3(A), 3(B) and 3(C). The substep 3(D) has 


only an aesthetic role, i.e., to have the general solution 


under the form: x: = fi(kı, ..., Kn-1), i= 1,n, fi being 
linear functions vith integer number coefficients. This 
"if possible" shows that substitutions are not always 


possible. But when they are we have to make all the 





possible substitutions. 

Note 1. The former algorithm is written under a form 
easy to introduce in the computer. 

Note 2. The former algorithm is more "rapid" than 
that of W. Sierpinski's 1, i.e., the general integer 
solution is reached after a smaller number of iterations 
(or, at least, the same) for any linear equation (2). In 
the first place, both aim at obtaining the coefficient + 1 


for at least one unknown variable. While Sierpinski 


started only by chance by decomposing the greatest 


coefficient in the module (writing it under the form of 


sum between a multiple of the following smaller 


coefficient (in the module) and the rest), in our 
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a 


algorithm this decomposition is not accidental but always 


seeks the smallest | x | and also chooses: the coefficients 


a; and a; for which this minimum is achieved. That is, 


test from the beginning the shortest way to the general 


integer solution. 


the shortest way; 


we 


Sierpinski does not attempt at finding 


the general integer solution of the equation and is not 


interested in how long it will be. However, when an 


algorithm is introduced in a computer, it is preferable 


that the computer time should be as short as possible. 


Example 1 


Let us solve in Z the equation: 17x - 7y + 10z = 


= -12. We apply the former algorithm. 


T h = 
2 r= 
3 |3| 
4 (1) 

y = 


= 1 


he knows that his way will take him to 
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Consider y: = ti h: = 2 


aş: = 3 p: = 2 





the other coefficients and variables remaining 


unchanged, go back to step 2. 








r = -1, i=l, j = 3 
l-1| =1 
x = —1(-3tz (-7tı) -12) = 3t, - 7ti - 12 
17-(-1) -1-17 
z = —1(17t> + (-7ti) ! + (-12)) = 
3 3 


-17t, + 42tı - 72 


(A) We substitute the values of x and z thus 
determined in the only relation (p) we 


have: 


(1) y tı + z = -17t2 + 43t: - 72 


(B) The substitution is not possible. 

(C) The substitution is not possible. 

(D) The general integer solution of the 
equation is: 


x = 3kı - 7k> + 12 


y = -17kı + 406572 
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z — —17kı + 42k> = T23 ki, kə 02 
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INTEGER NUMBER SOLUTIONS OF LINEAR SYSTEMS 
Definitions and Properties of the 


Integer Solution of a Linear System 


Let O aijXj = be, i — 1,m (1) 


A linear system with all its coefficients' integer 


numbers (the case with rational coefficients is reduced to 
the same). 


0 





Definition 1. Xj; = Xi, j = l,n is a particular 
0 Rə n 0 
integer solution of (1) if x 02, j= 1,n and O arx: = 
del 
bi, i m. Let the functions be f; BP 6 Z, j= dum 


Definition 2. Xj = f:? (kı, ..-kn), j = 1,n, is the 
general integer solution for (1) if: 





n E 
(a) Ó aif; (kı, əəə Kn) = bi, i = Lyn 

j=1 

of (ki, dE kn) 0 A; 

ğ ə. 
(b) For any Xj; = Xj, J = 1,n, particular 
0 0 
integer solution of (1), there is (ki, ..., 


0 
0:2 so that E (ki, es Sg kn) = X:, V; = lins 
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(In other words, the general solution is 


the solution that comprises all the other 


solutions.) 


Propriety 1 





A general solution of a linear system of m equations 
with n unknowns, r (A) =m <n is n - m times 


undetermined. 


Proof 

We assume by reductio ad absurdum that it is of order 
r, Pe ro «on mn, (the case r = 0, i.e., the Solution is 
particular, is trivial). It follows that the general 


solution is of the form: 
Xj > Ui1pi + ... + UirPr * Vi 
(Sı) 
Xn = UniP1 T gebe SE UnrPr + Vn, Uihr Vi 0 Z 


Pn = parameters 0 Z 





We prove that there are n m times undetermined 
solutions. The homogenous linear system (1), solved in r 


admits the solution: 


1 


1 
Dm Dn 
Xi ———— Xmeb o x >, 
D D 
Dm Dn 
Ra SS eds Vİ ea E Xn 
D D 
i É 
Let xi = xi, i = 1,n be a particular solution of the 
linear system (1). 
Considering 
Xmı = D ! Kmi, we get a solution 
1 1 0 
Xı = Dm+1Km+1 Test Dakn + Xı 
Xa =D l kn 
m m 0 
Xm — Dm+1Km+1 p Dakn + Xm 
0 
Xm+1 = D ! Km+1 2 Xm+1 
0 
Xn =D ! kan + x, kş = 
= parameters 0 Z 
which depends on the n - m independent parameters, 
system (1). 


for the 
Let the solution be n 


— m times undetermined: 
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Xı = Cııkı ün mi E Cin-mKn-m t dı 
(S2) 

Xn = Cniki dE əd Et Cnn-mEn-m + dn 

Cij, di 0 Z, Kj = parameters 0 Z 


(There are such solutions, we have proven it before.) 


Let the system be: 


aıXı + ... + aınXn = bı 

s aij, bi 02 

AmıXı + ... + amXn = Da Xi = unknowns 0 Z 
T: The case b; = 0, i = 1,m results in a homogenous 
linear system: 

aiiXi + ... + ain = 0, i= 1,m 
(S2) Y aii(CiiKi + ... + Cin-mKn-n + di) + ... aind(CaıKı + 


T: ege Cnn-mKn-m + da) = 0 
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0 = (a41C11 + ... + ainCnı)Kı + ... + (daiiCin-n + 


Fega T ƏinCnn-m) ! Kn-n + (aiidi “əə Əb 


+ ainda) , V k: 0 Z 


For ki = = ka-4 = 0 Y ajıdı + ... + aid, = O 

For kı = = kanı = kp = ... = Kin = 0 and ky = 
= 1 Y (acın + ... + aiCın) + (auda t... + 
+ aindan) = O Y âjıCın + ... + AinCnn = O 


V i= I,m, Vh = 1,n-m 








Cin | 
Vect. Vn = . , h = 1,n-m, are the particular 

Cnh solutions of the 
system. 
Vn, h = 1,n-m are also linearly independent because the 
solution is n - m times undetermined. {Vip vp Vana) cod 


is a linear variety that includes the solutions of the 


system obtained from (S2). Similarly, for (Sı) we deduce 


Uis 


that . SP S= l,r are particular solutions of 
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Uns 





the given system and are linearly independent because 





(Si) = r — times undetermined solution and 

Vi 
V = z =a solution of the given system. 

Vn 

The case (a). Ui, ..., Ur, œ = linearly dependent, 
it follows that (Ui, ..., Ur} is a free submodule of order 


r € n — mof solutions of the given system, then, it 
follovs that there are solutions that belong to (Vı, ..., 
Vn-m} + d and which do not belong to (Ui, ..., Ur}, a fact 
vhich contradicts the assumption that (Sı) is the general 
solution. 


The case (b). Ui, ..., Ur, V = linearly independent. 





(Ui, ..., Ur} + Vis a linear variety that comprises the 
solutions of the given system, which were obtained from 
(Si). It follows that the solution belongs to (Vı, ..., 
Vn-m} + d and does not belong to (Ui, ..., Ur} + V, a fact 
which is in contradiction with the assumption that (Si) is 


the general solution. 


II. When there is an i 0 1,m, b: .. 0 Y nonhomogeneous 


linear system 


(S2) 


aiiXi 


Y aii 


Tees 


(Ciiki + 


dAinXn 


+ Ain (Cnıkı + 


it follows that 


for 


for 


Y (aiiCi: + 


kı = 


kı = 


= Diy i= 1,m 


+ Cin-mEn-m + dı) + 


+ 


+ Cnn-nKn-m + dn = bi; 
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+ AinCni) Ki + ... + (@i1Cin-m + 

Tto... t+ AinCnn-m) Kim + (aıidı + + Aindn) 
= bi; 

= knm = Ü Y adı + ... + aindn = bi; 

= K3-1 = K 341 = ... = kam = 0 and kj = 1 Y 
Y (ai0Ci; + + AinCny + (aidi + ... + 
+ ain da) = bi it follows that 
t ... Tt AinCny = 0 
+ + aind bi, Vi= lm, VJ 1,n-m; 


49 


, j = 1,n-m are linearly independent 





because the solution (Sə) in 


n-m times is undetermined. 


di 


2?! Vs, j = 1,n-m andd = y are linearly independent. 


We assume that they are not linearly independent. It 


follows that 


d = sıVı + ... + Sn-mVn-m 

SiCim Pt ww T ÖnemCin-m 

SıCnı t ... + Sn-mCnn-m 
Irrespective of i = 1,m: 


bi = audi + 


ear əb ƏinÖn = dail (S1Cii ke dex Sn-nC1n-n) + 
11 
Te $e din (silin T ores FE SaumCuHnoan) = “agree 2. F 
+ ainCnı) Sı € ... + (@i1Cin-m + 


enə CE BinCnn-m) 


Sn-m 


Then, 


hypothesis 


bi = 0, 


irrespective of i= 1,m, 


contradicts the 


(that if there is an i10 1,m, bı... 0). It 


follovs that Vı, ..., Vn-n, d are linearly independent. 
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(Vi, ..., Vn-m} + d is a linear variety that contains 
the solutions of the nonhomogeneous system, solutions 
obtained from (S2). Similarly, from (Sı) it follows that 
{Gi, ..., Gr} + V is a linear variety containing the 
solutions of the nonhomogeneous system, obtained from 
(Sı). 

n-m > r follows that there are solutions of the 
system that belong to (Vı, ..., Vnm} + d and which do not 
belong to (Gi, ..., Gr} + V (it contradicts the fact that 
(Si) is the general solution). Then, it results that the 
general solution depends on the n-m independent 
parameters. 

Theorem 1. The general solution of a nonhomogeneous 
linear system is equal to the general solution of an 
associated linear system plus a particular solution of the 


nonhomogeneous system. 


Proof 


Let the homogeneous linear solution: 


ƏıXı + ... + aınXn = 0 
, (AX = 0) 


ƏmıXı + ... + AmXn = Ü 


with the general solutions: 


Xı = Cııkı Web E Cin-mKn-m + di 
Xn = Cniki To em SE Cnn-mEn-m F dn 
0 
and Xı = Xı 
. 0 
Xn = Xn 


a particular solution of the nonhomogeneous linear 


0 
çi Xı = CııKı + ... + Cin-mKn-m + d + Xi 


0 
XQ — Ontki P ə CBOC€anKsew CE dua Geox 


is a solution of the nonhomogeneous linear system. 


We have written 


Alı ... dın Xı bı o 
A= : xD E işe Əə ə , OF. 
Ami d inn Xn Da o 


system 


52 


(vector of dimension m), 


ki Cil ... Cin-m dı 
Xı 
K = , C= , do 5 
; 
.0 
Kass Cubo xum dn 
Xn 
AX = A(Ck + d + x) = A(Ck + d) + AX =b+02=b 


We will prove that irrespective of x1 = yi 


. 0 
Xn = Yn 


there is a particular solution of the nonhomogeneous 
system 


kin = Kim 0 Z with the propriety: 


0 0 0 0 
Xi  Cııkı + ... + CinKn-m + di + Xı = yı 


. 0 0 0 0 
Xn = Cniki mu 1 Cnn-mEn-m + di + Xn = Yn 
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We write y° = yı 
-0 
Yn 
Ra 2 . e 
We demonstrate that those k; 0 Z, j = 1,n-m are those 
for 
0 
which A(CX? + d) = 0 (there are such X, 0 Z because 
Xi = 0 
Xn = 0 


is a particular solution of the homogenous linear system 


and X = CK + d is a general solution of the 
nonhomogeneous 


linear system) A(CE” + d + X? - Yİ) = A(CK? + d) + AX? - 
- A? =0+b-b=0. 


Propriety 2. The general solution of a homogenous 





linear system can be written under the form: 


Xı € Ci1ki + ... + Cin-mKn-m 


Xn = Cnıkı + ... + Cnn-mEn-m , 
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k: = a parameter belonging to Z (with di = d, = = 
= dn = 0) 
Proof 
(SG) = general solution. It results that (SG) is 


(n-m) times undetermined. 


Let (SG) of the form be 


Xj = Cııbı E eres F Cin-mPn-m Li dı 
Xn = Cnipi F a T Cnn-mPn-m iLE dn 
with not all d; = 0; we demonstrate that it can be written 
under the form (2); the system admits the trivial solution 
Xı = 0 0 Z 
; 
Xv > 002 


it results that there are p; 0 Z, j = 1,n-m 


0 0 
Xj = Cııbı "E əə, Ub Cin-mPn-m + dı € 0 
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(4) 
. 0 0 
Xn = Cnipi "ou ə tt Cnn-mPn-m + da = 0 


p; OZAYp; 0 2 


which means that they do not make any restrictions. 


It results that 


0 0 


Xı = CııKı + ... + Cin-mKn-m + (Ciipi Fo... Cc Cin-mPn-m + di) 
. 0 0 
Xn = Cniki Bole E C'un-mKn-m + (Capi F owe OF Cnn-mPn-m F dn) 
0 0 əz 
But Cpıpı + ... + Chn-mÖn-n + di = O, h = 1,n (from (4)). 





Then the general solution is of the form: 


žr daa by b “Gina Ee 


eee + Cnn-mEn-m 


kj; = parameters 0 Z, j = 1,n-m, it results that dı = dz = 


Theorem 2. Let the homogenous linear system be: 


ek +... + AainXn = 0 
— +o... + amXn = 0, r(A) = m; 
(ahı, ..., Ann) = 1, h= Lu and the general solution 
2 = Ciiki + ... + Cin-mKn-m 
3 = Caıkı 4 22. CF Con-mkn-m 
then (anı, ..., dhi-i, Anitiy «eer ann) | (Grecs inm) 
irrespective of h - Ton and i = 26 


Proof 


Let some arbitrary be h 0 1,m and some arbitrary i 0 


eee + aàni-i1Xi-1 + ƏhasıXauı + ... + AnnXn = AniXi- 
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Because (ahı, ..., dhi-i, A@nitiy «++, Ann) | ani it results that 
d = (ahı, ^, Ani-ly Öhiil: 55-57, Ann) b” irrespective of 


the value of xi in the vector of particular solutions; for 


k2 = k3 = ... = knm = Ü and ki = 1 we get the particular 
solution: 

Xj = Cii 

Xi = Ci Y d| cis and so on 

Xn = Cni 
for Xı = k2 = ... = kimi = 0 and krm = 1 the following 


particular solution results: 


X1 = Cin-m 
Xi = Cin-m 
Xn = Cnn-m 
it results that dccus hence, d| cis, j = 1,n-m Y 


Y d| (City ..... Cin-m) . 
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Theorem 3. 
Xi = Cııkı + ... + Cin-mKn-m ci, 0 Z being given 
If 


Xn = Cniki + ... + Cnn-mKn-m, Kj = parameters 0 Z 


is the general solution of the homogenous linear system 


ƏıXı + ... + AinXn = O m € n, 

amıXı + ... + amön = Ü r(A) =m; 
then (cı), ..., Cn) = 1, V j = 1,n-m. 
Proof 


We assume, by reductio ad apsurdum, that there is 


20 0 1,n-m: (Co 4 ET Cni, ) =d.. 1, we consider the 





maximal co-divisor > 0; we reduce the case when the 


maximal co-divisor is -d to the case when it is equal to d 
(nonrestrictive hypothesis); then the general solution can 





be written under the form: 


Xı = Ciiki + ... + C1j dk; FP oue c Cini Kaen 
0 0 
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(5) ; 
Xn = Cakı + ... + Cnj dk; Fo ve oC pnt ka 
0 
where d = (ciy / pepe 00 Cag sd oig cand (Cii ; Ar. 
e... Cnj ) ^-^ 1. 
0 


We prove that 


is a particular solution of the homogenous linear system. 


We write 


Cri sə Cij 2 Sl xu Cin-m kı 
C= . > . k = . ; 
Cni ... Crj d ... Can-m kj 
0 0 
Kn-n 
x =c ! k the general solution; 
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We know AX = 0 Y A(CK) = 0, A= 


We assume that the principal variants are xi, wa. Xm (LİE 
not we have to renumber). It follows that Xmı, ... Xn iS 
the secondary variant. 


For kı = e.s. = Kj = ki = e.s. = Kn-m = 0 and ks = 1 


we get a particular solution of the system 


Xı = Cij d C1j d C1j 
. 0 
Y 0-A -dlA Y 
Xn = Cnj d Ora. d Grae 
0 
C13 Xı = Cij 
. 0 
Y A = 0 Y 
Cnj Xn = Cnj 


is the particular solution of the system. 
We demonstrate that this particular solution cannot 


be obtained by 


Xı = Ciiki + ... C1j dk; Toe e E CinmKi-m =” CHI 


(6) 
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Xy — eyra, A seeps: FF ə dk; İmar GH Chaska = Cnj 
Xmt1 = Cmıkı + ... + Cnadk; que şir Ys Cms1,n-mKn-m — Cnáj 
(7) 
Xn = Oniki bı OF Caj dk; Poanse Fo naemKnsm cm o 
Cm+1,1 e+. Cm+1 4 ... Cm+in-m 
0 
Chi .. Cn3 “ə. Chien 1 
Y ky = = - 6 Z (because 
d..1) 
0 i d 
Cmt+1,1 e+. Cm+1, 5 GQ) x Giri, ns 
Cni «- Ong «Gd. ws Cn,n-m 


It is important to point out the fact that those kj = 
ki, 


j = 1,n-m that satisfy system (7) also satisfy system 
(6), 


because, otherwise (6) would not satisfy the definition 
of the solution of a linear system of equations (i.e., 
considering system (7) the hypothesis was not 


restrictive). From X; 0 Z it follows that (6) is not 


the 
0 


general solution of the homogenous linear system 


contrary to the hypothesis); 


irrespective of j = 


then (cı: Cani) = 





Propriety 3. 





aixi + ƏınXn 


ƏmıXı + 
aij, bi 0 


Z, r(A) - 


Solved in R, 
£1 (Xm, n 
Xn = 


£n (Xm, .... 


main variants vhere f: 


+ amin = 


m € n, 


Let the linear system be 


= bi 


Dn 


xj; = unknowns 0 Z 


we get 


Xn) 


eK. daş Xn are the 


Xn) 


are linear functions of the form: 


gi —-— 





d: i 
Cm+1Xm+1 + . . +CnXntei 
fi = 
di 
j = 1,n-m 
Ci 
If —— 0 Z irrespective of i = 


where Cmj, di, e1 0 2; i = 


1,m then the linear system 
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di 


admits integer solution. 


Proof 


For 1 < i € m, Xi 0.2, then fj 0 Z. Let: 


Xmtı = Umtikm+1 
Xn = unkn 
Solution 
1 1 Cı 
Xı VmsıKmt1 + + Vaka T Pod 
di 
m m Em 
Xm = Vari Emea E e. E VnKn + Pn 
Am 


where Umı is the maximal co-divisor of the denominators 
of 


i 
Cm+j — MN n 


the fractions --- , i = 1,m, j = 1,n-m calculated after 
di 


their simplification when they were irreducible. 





Vig: “SS O 2; this is a solution n-m times 





undetermined (depends on n-m independent parameters: 


Kmeny ---., Kn) but is not a general solution. 
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Propriety 4. 





if 


Under the conditions of propriety 3, 


io io 


there is an io O 1,m: fi = UmtiXme1 + ... + Un Xn + —— 
0 
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ó Z then the system 


If there 


Ci 
x 0 
yer, e o SS 
with Um; O Z, j = 1,n-m and -— 6 
does 
di 
0 
not admit integer solution. 
Proof 
© Xm; ..., Xn in Z it results in Xi, ö 2. 
Theorem 4. Let the linear system be 
aııXı + ... + aınXn = bı 
AnıXı + ... + amXn = Dn 
diş, b: 0 Z, Xj = unknowns 0 Z, r(A) =m < 
are 


indices 1 < i, < € İn € n, 


h = 1,m with the propriety: 


mn 


in 0 (1, 


and 


2, 


eg. 


n], 
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E = bi aii e aii 
] 2 m 
11 
is divided by 
Dm Ami Ami 
2 m 
ali aii bi 
a x. ı m-1 
x m . . 
İn : i is divided by 
Ami Ami bn 


then the system admits integer number solutions. 


Proof 


We use propriety 3 


di = u^ L = 1,m; ei = Pie , h = 1,m 
h an 
Note 1. Conversely, it is not true. 
Consequence 1. Any homogenous linear system admits 





integer number solutions (beside the trivial one); r(A) 


Proof 


a 


ex = 0 ! ?, irrespective of h= 1,m. 


Consequence 2. If ê= + 1, it follows that the 





linear system admits integer number solutions. 


Proof 


a I (+1), irrespective of h = 1,m; 
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FIVE INTEGER NUMBER ALGORITHMS TO 


SOLVE LINEAR SYSTEMS 


This chapter further extends the results obtained in 
4 and 5 (from linear equations to linear systems). Each 
algorithm is strictly demonstrated and then an example is 
given. 

Five integer number algorithms to solve linear 


systems are further given. 


Algorithm 1 (method of substitution) 





(Although simple, this algorithm requires complex 
calculus but is, nevertheless, advantageous in introducing 
it in the computer). 

Some integer number equations are initially solved 
(which is usually simpler) by means of one of the 


algorithms 4 or 5. (If there is an equation of the system 





which does not admit integer systems, then the integer 
system does not admit integer systems. Stop.) the 
general integer solution of the equation will depend on n- 


1 integer number parameters (see 5): 


( 
(pı) Xi = fi (ki, ww Kn-ı ), s]. m l,n where all 
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(1) 
the functions f: are linear with integer number 
1 


coefficients. 
This general integer number system (pı) is introduced 


in the other m-1 equations of the system. We get a new 





system of m-1 equations with n-1 unknown variables: 


ki , ii = 1,n-1, which is also to be solved in integer 
numbers. The procedure is similar. Solving a new 
equation, we obtain its general integer solution: 


(1) (2) (2) (2) f 
(p2) ki = fi (Ki , ..., Kn-2 ), i2 = 1,n-l 


where all the functions R are linear, with integer 
number coefficients. (If, along this algorithm we come 
across an equation which does not admit integer solutions, 
then, the initial system does not admit integer solution. 
Stop.) 

In the case that all the solved equations admitted 
integer systems at step (j), 1 < j < r, (r being of the 


same rank as the matrix associated to the system) then: 


(3-1) (3) (3) (3) 
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(ps) ki = Eu (kı , , Kn-j ), Y = l,n-jJtl, 
(3) 
fi are linear functions vith integer number 
2 
coefficients. 


Finally, after r steps, and if all the solved 
equations admitted integer solutions, we get to the 
integer solution of one equation with n-r+1 unknown 
variables. 

The system will accept integer solutions if and only 


in this last equation will have integer solutions. If it 





does, let the general integer solution of it be: 


(r) 
where all fi are linear functions with integer number 
T 


coefficients. 


Now the reverse way follows. 


(sty Ue 


We introduce the values of ki y i. = 1,n-r+1 
x «Garner 
at step (pr) in the values of 5. g? i-1 = 1,n-r+2 from 
step (Drı). 


It follovs: 
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(r-2) (r-1) (r) (r) (r) (r) (r) (r 


) 
ki = fi (fi (ki esə Kezer eş Enes (ki peeerKn-r) ) 
m r-1 r-1 
(r-1) (r) ©. tg TIS NIS 
= Gi > (kı , , Kn-r ), 1r-1 = l,n-r-1 
(=I) 
from which it follows that gi are linear functions with 


integer number coefficients. 
Then follow those from (pr-2) in (pr-), and so on, 


until we introduce the values obtained at step (p2) in 


those from the step (pi). It will follow: 
(1) (r) (r) 
Xi, = gi (ki , n Kicz notation Ji, (ki, . Kn-r) , 
is işn 


with all gi T most obviously, linear functions with 
integer number coefficients (the notation was made for 
simplicity and aesthetical aspects). This is, thus, the 
general integer solution, of the initial system. 


The correctness of algorithm 1. The algorithm is 





finite because it has r steps on the first way and r-1 
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steps on the reverse. (r € + 4). Obviously, if one 
equation of one system does not accept (integer number) 
solutions then the system does not have either. 

Writing S for the initial system and S: the system 
resulted from step (pj), 1 € j € r-2, it follows that 
passing from (p;) to (pa) ve pass from a system S: to a 


system S311 equivalent from the viewpoint of the integer 
(3-1) Be ” ea ERN e 

number solution, i.e., Ki = t; , iş € 1,n-j*1 which is 
3 2 


a particular integer solution of the system S: if and only 


if kk = hi, , iss. = 1,n-) is a particular integer 


0 (331) 0 0 
solution of the system Ss: where hi = fi (ti, 6.8, tne 


jn) , . . 
TEL j*1 


ini = 1,n-j. Hence, their general integer solutions 





also equivalent (considering these substitutions). So 
that, in the end, the solving of the initial system S is 


equivalent with the solving of the equation (of the 
system 


consisting of one equation) Sr-1 with integer number 


coefficients. It follows that the system S admits 
integer 


number solution if and only if all the systems S; admit 
integer number solution, 1 € j € r-1. 


Example 1. By means of algorhythmus 1, let us 





calculate the integer number solution of the system: 


5x — 7y — 2z + öv = 6 
(S) 
9 -4x + öy — 3z + ll = 0 


Solution: We solve the first integer number equation. 


We 


obtain the general integer solution (see [4] or [5]): 


X = ti + 2to 


Mec En 

(p1) 
z = -tı + 5t2 + 3t3 - 3 
w = ts , 


where tı, to, ts 02. 


Substituting in the second, ve get the system: 


(Sı) 5ti — 23t> + 2t3 + 9 = 0 
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Solving this integer equation ve obtain its general 


integer solution: 


tı = ki 
(p2) to = kı + 2k» + 1 


t3 = 9ki + 23k, + 7 , 
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where ki, Kə 02. 


The reverse way. Substituting (pz?) in (pi) we 





obtain: 
x = 3ki + Ak, + 2 


ki 


K 
I 


z = 31k; + 79k. + 23 


w = 9k + 23k: + 7 , 


where ki, k2 0 Z which is the general integer solution 
of 


the initial system (S). Stop. 


Algorithm 2 





Input 


A linear system (1) without all ai; = 0. 


Output 
We decide on the possibility of an integer solution 


of this system. If it is possible, we obtain its 
general 


integer solution. 


Method 


1. t-1,h-1,p-1 


unknown 


the 


occurs, 


identity, 


equations 


there 


76 


(A) Divide each equation by the maximal co- 


divisor of the coefficients of the 


variables. If you do not get an integer 


quotient for at least one equation, then 
system does not admit integer solutions. 


Stop. 

(B) If there is an inequality in the system, 
then the system does not admit integer 
solutions. Stop. 


(C) If the repetition of more equations 
keep one and if an equation is an 


remove it from the system. 
If there is (io, jo) so that lass. | = 1, then 
obtain the value of the variable Xi, from the 
equation io; relation (Ti). Substitute this 


relation (where possible) in the other 


of the system and in the relations (Tt-1), (Hn) 


and (Pp) for all i, h and p. Consider t: = t 
remove equation (io) from the system. If 


is not such a pair, go on to step 5. 


4. 
unknovn 


and 


minimum 


variants, 


6. 
all 
relations 

7. 
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Does the system (left) have at least one 
variable? If it does, consider the new date 


go on to step 2. If it does not, write the 
general integer solution of the system 
substituting ki; ko, ... for all the variables 
from the right term of each expression which 
gives the value of the unknowns of the initial 
system. Stop. 


Calculate a = min tlel, aii, / r (mod ais), 
Ty ji: İz 





o< |r| < |as; |, and determine the indices 


i, jı; İz as well as the r for which this 


can be calculated. (If there are more 





choose one, arbitrarily.) 


Write: xj = tı ——————— Xi; , relation (Hn). 


Substitute this relation (where possible) in 


the equations of the system and in the 


(Te), (Hn) and (Pp) for all t, h and p. 


(A) If a . 1, consider x; : = tn, h: =h * 1 








and 
2 
go on to step 2. 
(B) If a = 1, then obtain the value of x; 
from 
1 
from the equation (i); relation (Pj). 
Substitute this relation (where possible) in 
the 
other equations of the system and in the 
relations (Tt), (Hn) and (Pp-1) for all t, h 
and p. 
Remove the equation (i) from the system. 
Consider h: = h + 1, p: =p + 1 and go back to 
step 4. 
The correctness of algorithm 2. Let the system (1) 
be. 
Lemma 1. We consider the algorithm at step 5. 
Also, 
let M= {|r|l, aij / r (mod aij v 0 < |e] < lau l; Ly 
Jı, 1 2 2 
> = ly 25 94. ee) Then Mw 
Proof 
Obviously, M is finite and Md N*. Then, M has a 
minimum if and only if M... We suppose, conversely, 
that M = . Then aj 27 0 (mod ai; ), V i, ji, j2. It 


follows conversely as well that aij, / 0 (mod ais ); V i, 


ji; ja. That is | ais | = | ais |, V Ly ji; ja. We 


consider 
1 2 


a io arbitrary but fixed. It is clear that (aii; NOE 
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.; din) — aij = 0, V j (because the algorithm has 
passed 
0 0 
through the substeps 2(B) and 2(C)). But, as it has 
also 
passed through step 3, it follows that lİa::l . 1, V j, 
but 


0 


as it previously passed through step 2(A), it vould 
result 


that la: s) = 1, V j. This contradiction shows: that the 


assumption is false. 


Lemma 2. Let ai, / r (mod aig). Substitute X 

Sd AS 

x o Xi, in the system (A) obtaining the system 
ur 

" DE 
(B). Then, x; = x3, j = 1,n is the particular integer 
0 

solution of (A) if and only if x; = x, j = Jo and th = 

m = 
0 
= 300 hs SSS d is the particular integer solution of 
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Lemma 3. Let a; . 1 and az be obtained at step 5. 
Then 0 € az € ai. 


Proof 


It is sufficient to show that a; € |aij|, V i,j 
because 


in order to get az step 6 is obligatory, when the 
coefficients of the nev system are calculated, a: being 
equal to a coefficient from the nev system (equality of 


modules), the coefficient on (io2ı). 


Let ai, with the propriety lai j | < Bi. Hence, aı 
> 00 00 
> la: 5| = min Class Let aij with laij | > | ais | 
there 

0 0 0 s 0 s m 





is such an element because İa:: | is the minimum of the 
0m 


coefficients in the module and not all |ai j|, j = 1,n 


are 
0 
equal (conversely, it would result that (ai 5, ..., aix) 


0 0 


— aij, V j 0 1,n; the algorithm passed through substep 
2(A) 
0 


has simplified each equation by the maximal co-divisor 
of 


its coefficients; hence, it would follow that la: s) = 1, 





V j = 1,n which, again, cannot be real because 


the algorithm has also passed through step 3). Of the 
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coefficients ai; we choose one with the propriety aij 

Om Os 

0 
581 4 there is such an element (contrary, it would 
result 
(aii, ... ain) = laij |, but the algorithm has also 
passed 
0 0 Om 
through step 2(A) and it would mean that |a; | = 1 
which 
Om 


contradicts step 3 through which the algorithm has also 
passed). 


Considering do = ai; / ai; O Z and r = a: - 
0 s Om Os 
z 0 aS 0 


— qo ai; 0 Z, we have ai j / ro (mod ai; ) and Ü € 
0m O's om 
0 


< İrol € lar: | € lar: | € ài. We have, thus, obtained 


an 
Om 00 


ro with İrol € ai, which is in contradiction with the 
very 


definition of a;. Thus, a: € İasİ, V i,j. 


Lemma 4. Algorithm 2 is finite. 


Proof 

The functioning of the algorithm is meant to 
transform a linear system of m equations and n unknowns 
into one of mx nı, with m; € m, nı € n and, thus, 
successively into a final linear equation with n - r + 1 


unknowns (where r is the rank of the associated matrix). 
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This equation is solved by means of the same algorithm 


(which works as [5]). The general integer solution of the 
system will depend on the n - 1 integer number independent 
parameters (see 161--similar proprieties can be 


established also for the general integer solution of the 
linear system). The reduction of equations occurs at 
steps 2, 3 and substep 7(B). Steps 2 and 3 are obvious 


and, hence, trivial; they can reduce the equations of the 





system (or even put an end to it) but only under 
particular conditions. The most important case finds its 


solution at step 7(B), which always reduces one equation 





of the system. As the number of equations is finite, we 
come to solve a single integer number equation. We also 
have to show that the transfer from one system mix nı to 


another mi+ixNi+1 is made in a finite interval of time: by 





steps 5 and 6 permanent substitution of variables are made 
until we get to a = 1 (ve get to a = 1 because, according 
to lemma 3, all a - s are positive integer numbers and 
form a strictly decreasing rov). 


Theorem of correctness. Algorithm 2 correctly 





calculates the general integer solution of the linear 





system. 


Proof 
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Algorithm 2 is finite according to lemma 4. Steps 2 
and 3 are obvious (see also [4], [5]). Their part is to 
simplify calculus as much as possible. Step 4 tests the 


finality of the algorithm; the substitution with the 
parameters ki, kz, ... has systemization and aesthetic 
reasons. The variables t, h, p are counter variables 
(started at step 1) and they are meant to count the 
relations of the type T, H, P (numbering required by the 
substitutions at steps 3, 6 and substep 7(B); h also 
counts the new (auxiliary) variables introduced in the 
moment of decomposition of the system. The substitution 
from step 6 does not affect the general integer solution 
of the system (it follows from lemma 2). Lemma 1 shows 
that at step 5 there is always a, because . Md N*. 

The algorithm performs the transformation of a system 
mixni into another, ai4axnia4, equivalent to it, preserving 
the general solution (taking into account, hovever, the 


substitutions) (see also lemma 2). 





Example 2. Calculate the integer solution of: 


— 12x — Ty + 92 12 


I 
o 


— 5y + öz + 10w 
Oz + Ow = O 


15x + 21z + 69w 


ll 
Cə 


Solution 
We apply algorithm 2 (we purposely looked for an 


example to be passed through all the steps of this 


algorithm: 
Ii. t s1, hzl, p < 1 
2. (A) The fourth equation becomes: 5x + 7z + 
23w = 1 
(B) teme 
(C) Equation 3 is removed. 
Bis No; go on to step 5. 
5. a = 2 and i = 1, ji = 2, j2 = 3 and r = 2. 
6. z = tı + y, the relation (Hi). Substituting it 


in the system: 


— 12x + 2y + 9ti = 12 
3y + 8ti + 10w = 0 
5x + 7y + 7tı + 23w = 1 
he a . 1; consider z: = ti, h: = 2 and go back to 
step 2. 
2. SS 
3 No Step 5. 


Bu a = 1 and i = 2, ji = 4, iz = 2 and r “1. 
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y = t, — 3w, the relation (Ho). Substituting in 


the system: 


— 12x + 2tə + 9t, - öv = 12 
3t> + 8tı + w= 0 
5x + 7t, + 7ti + 2w = il: 


Substituting it in relation to (Hı), ve get: 


z tı + te — 3w, relation (Hı)". 


mi — 3te - 8ti, relation (Pi). 
Substituting it in the system, we get: 
— 12x + 20t2 + 57ti = 12 


5x + te - 9ti = 1 


Substituting it in the other relations, we get: 


z = 10t2 + 25ti , (Hi) ''; 

y = 10t> + 24ti , (Hə) "1, 

h: = 3, p: = 2 and go back to step 4. 
Yes 


te = 1 — 5x + 9tı, relation (Tı). 
Substituting it (where possible) we get: 


{- 112x + 237ti = -8 (the new system); 


z = 10 - 50x + 115ti, (Hi) ''! 
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y = 10 — 50x + 114t;, (H2) '' 

w= -3 + 15x - 35tı, (Pı) " 

Consider t: = 2; go on to step 4. 

Yes. Go back to step 2. (From now on algorithm 


2 works Similarly with that from [5], with the 
only difference that the substitutions must also 
be made in the relations obtained up to this 
point). 


No. Go on to step 5. 


a = 13 (one three) and i € 1, Jı = 2, jo = 1 and 
r= 13. 
x = t3 + 2ti, relation (H3). 


After substitution we get: 


-112t: + 13ti = -8 (the system) 
z = 10 = 50t; + 156: (H)™;}; 


y = 10 - 50t3 + 14ti , (H2) "r 


w = —3 + 15t:- Sti y (Pi) '! 
tos. Sİ. to tı , (Ti)! 
x: = t3, h: = 4 and go on to step 2. 


No, go on to step 5. 


54 a=5 and i = 1, jı = 1, j2 = 2 and r ^ 5 


6. tı = ta + 9t3, relation (Hi). 
Substituting it, we get: tə + 13t, = -8 (the 
system). 


z = 10+ 85t4 + 15tu , (Hi)” ; 


y = 10 + 76t3 + 14t, , (H) 25 


x = Tote A, “tan. p- (Ha) ou 
w = -3 — 30t3 - Sta , (Pa) t Phy 
t2 = 1 - 14t; - ta r Alak 7 
Tas ti: = ta, h: = 5 and go back to step 2. 
2. m 
3 No; step 5. 
Bix a= 2and i € 1, jı = 2, j> = 1 and r / -2. 
6. t3 = ts — 3t,, relation (Hs). After 
substitution, 
we get: 
5t: — 2t, = -8 (the system) 
z = 10 + 85ts - 240t, , iy Vg 
y = 10 + 76t5 - 214t, , (H2) ; 
x - 19t5 — 55ta, (Hs)”” ş 
w = -3 - 30ts + 85t4 , (Pa) 3 
t2 = -l — 14ts + 41t. , (Ti) '''; 
ti = 9ts + 26t, , (Ha) ' ; 





Ta tai = te, 
2. E 
3 No; step 5. 
Əs a = 1 and i = 1, jı = 
6. tu = te + 2t:, relation 
substitution, 
we get: 
besa 2 EE -8 
Zz = 10 - 395t: - 240t: 
y = 10 - 392t: - 214t: 
x = — TES — 
w = -3 + 140ts + 
t2 = + 68t: + 
tı — = 43t: - 
ta = = Sts) = 
JT ts = 2t, — 8, relation 


in 


h: = 6 and go back to step 2. 


the system, we get: 


S5t. 
85t6 
Alte 
26t6 


Bte 


0 


(He). 


ja: 


(P2). 


Ox 


r= 1. 


After 


(the system) 
(H) 3 
(H) 4 
(Ha) t"! 

(Pi)" 3 
(T) q 


(Ha) '' ; 


(Hs) ' ; 


Substituting it 


Substituting it in the other relations, it 


follows: 


N 
I 


K 
I 


—-1030t6 + 3170 


918t. + 2826 
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x = — 237te + 
w = 96btç, = 
t2 = 177t: - 
ta = 112t: + 
importance 
t3 = 13te + 
ta = Ste — 
Consider h: = 7, 
te 0 Z 
4. No. The general 
system 
is: 
x = —-237k,ı + 728 
y = —918ki + 2826 
z = 1030k: + 3170 
w= 365ki - 1123 


728 
1123 
543 


344 
relations of no 


p: = 3 and go back to step 4. 


integer solution of the 


where ki is an integer number parameter. 


Stop. 


Algorithm 3 





Input 


A linear system (1). 


Output 
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90 
We decide on the possibility of an integer solution 
of 
this system. If it is possible, we obtain its general 


integer solution. 


Method 
T. Solve the system in R". If it does not have 
Solutions in R^, it does not have solutions in 
2^ 
either. Stop. 
2 f=1, t=1, h=1, g=1 
SG Write the value of each main variable x; under 
the form: 
(Ee,i)it Xi = Ó qij X; t qi + (Ó rx: unb 
J J 
With all gis; dir Tu, Ifa; an Z so that ali 
İml < Jal, & 200, İril € |&| (where all x; 
of 
the right term are integer number variables: 
either of the secondary variables of the 
system 
or other new variables introduced with the 
algorhythmus). For all i, we write ris yem 
4. (Fe,i)a: Ö eke = fi yYr,i + ri = 0 where (Ye,i)i 


J 


of 


the 


expression 


the 


and 
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are auxiliary integer number variables. We 
remove all the equations (Fr,ı) which are 
identities. 

Does at least one equation (Fr,ı) exist? If it 


does not, write the general integer solution 
the system substituting ki, ko, ... for all 
variables from the right term of each 


representing the value of the initial unknowns 

of the system. Stop. 

(A) Divide each equation (Fr,ı) by the maximal 
co-divisor of the coefficients of their 
unknowns. If the quotient is not an 


integer number for at least one iç then 
system does not admit integer solutions. 


Stop. 

(B) Simplify--as in m--all the fractions from 
the relations (Er,i)i. 

Does it exist Ti jig with the module 1? 

If it does not, go on to step 8. 


If it does, find the value of x: from the 


equation (Fe ); write (Tt) for this relation 


t, 


Consider 


variants, 


10. 


0 


substitute it (where it is possible) in the 


relations (Eşi), (rt-1), (Hn), (Gg) for all i, 


h and g. Remove the equation (Fs, ). 





0 


f: = f + 1, t: =t + 1 and go back to step 3. 


Calculate a = min tlel, Li, / (mod ri), 
Ly ji: İz 





0 |r| < [ri 1) and determine the indices in, 


ji, İz as well as the r for which this minimum 





can be obtained. (When there are more 


choose only one). 


(A) Write x; = Zh — —————— Xj; , where zy is 


new integer variable; relation (Hi). 


(B) Substitute the letter (where possible) in 





the relations (Egi), (Frei), (Te), (Hn), 
(Gg) for all i, t, h and g. 


(C) Consider h: = h + 1. 


(A) If a. 1, go back to step 4. 


(B) If a = 1, calculate the value of the 


variable x; from the equation (Fs); 
1 
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relation (Gg). Substitute it (where 


possible) in the relations (Esi), (Tt), 
(Hn) r 


(Gg-1) for all i, t, h and g. Remove the 
equation (Fgi). Consider g: = g + 1, 


f: = f + 1 and go back to step 3. 


The correctness of algorithm 3 





nə ' 
Lemma 5. Let i be fixed. Then ( 0 rix: + vi) /5: 
dn: 


(with all ru, rip “iç ni, m? being integers, ni < ns, ñi a: 


0 


and all x: being integer variables) can have integer 
values 


if and only if (Lin ep Tin, / Bef us 


Proof 


The fraction from the enunciation can have integer 


nə 
values if and only if there is a z 0 Z so that ( Ö 
rijXj + 
j=n: 
no ' 
+ 14)/4 =z] Ó rix? - % z * ri = 0 which is a linear 
den: 
equation. This equation admits integer solution 


] (Xin , wswp Piny a) ses. 
1 2 
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Lemma 6. The algorithm is finite. It is true. The 
algorithm can stop at steps 1, 5 or substep 6(A). (It 
rarely happens to stop at step 1). An equation after 


another are gradually eliminated at step 4 and especially 
7 and 10(B) (Fgi)--the number of equations is finite. If 
at steps 4 and 7 the elimination of equation may occur 
only in special cases, elimination of equations at the 


substep 10(B) is always true because, through steps 8 and 


9 we get to a = 1 (see [5]) or even lemma 4 (from the 
correctness of algorithm 2). Hence, the algorithm is 
finite. 

Theorem of Correctness. The algorithm 3 correctly 





calculates the general integer solution of the system (1). 





Proof 
The algorithm is finite according to lemma 6. It is 


obvious that if the system does not have solutions in R" 





it does not have in Z" either, because Z" d R" (step 1). 
The variables f, t, h, g are counter variables and are 
meant to number the relations of the type E, F, t, H and 
G, respectively. They are used to distinguish betveen the 
relations and make the necessary substitutions (step 2). 
Step 3 is obvious. All the coefficients of the unknovns 


being integers, each main variable x: vill be vritten: 
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Xy = (Ó Cika caer y 
j 


which can assume the form and conditions required in this 
step. Step 4 is obtained from 3 by writing each fraction 
equal to an integer variable yr,: (this being xi 0 2). 

Step 5 is very close to the end. If there is no fraction 
among all (Er,ı) it means that all the main variables xi 
already have values in Z, vhile the secondary variables of 
the system can be arbitrary in Z, or can be obtained from 
the relations T, H or G (but these have only integer 
expressions because of their definition and because only 
integer substitutions are made). The second assertion of 
this step is meant to systematize the parameters and 
renumber, it could be left out but aesthetic reasons 
dictate its presence. According to lemma 5 the step 6(A) 
is correct. (If a fraction depending on certain 
parameters (integer variables) cannot have values in Z, 
then the main variable vhich has in the value of its 
expression such a friction cannot have values in Z either, 
hence, the system does not admit integer systems). This 
6(A) also has a simplifying role. The correctness of step 


7, trivial as it is, also results from [4], and the steps 
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8-10 from [5] or even from algorhythmus 2 (lemma 4). 

The initial system is equivalent to the "system" from 
step 3 (in fact, (Eri) as well, can be considered a 
system). So, the general integer solution is preserved 
(the changes of variables do not prejudice it (see [4], 
[5], and also lemma 2 from the correctness of algorithm 
2)). From a system mxn;i we form an equivalent system 
MişıXNişı With mu: € m; and nu: € ni. This algorithm works 
similarly to algorithm 2. 

Example 3. Employing algorithm 3, find an integer 


solution of the following system: 


3x: + 4x5 + 22x, — 8x5 = 25 
6X1 + + 46x, — 12x: = 2 
4X» + 3X3 = X4 + 9x5 = 26 
Solution 
1. Common solving in R”, it follows: 
23X4-6X5—1 
xt SSS a 
— =3 
X4T2Xx5424 
KS o SS 
4 
11x5+2 
X30— mi 
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2x4-1 
Xi = -TX4t2X5 + ———— 
- (Er,ı) 


5. 


o 


d 
relation 
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X4T3X5 
Xx Lp exces (E1,2) 
4 
X52 
X3 —4Xs a (E1,3) 
3 
2X4 + 3yii — 1 - 0 (Fı,ı) 
X4 + 2X5 — 4yı2 - 0 (F1,2) 
Xs — 3yis +2=0 (Fı,3) 
Yes 
Yes [rss] = 1. Then X5 = 3yis = 2, the 
(Tı) Substituting it in the others, ve get: 
2x4-1 
Xı = —x4t 6yı3>-4 t pi ——--— 
-3 
X4t6y13—-4 
ax Ge NOUO Ue (E1,2) 
4 
3yi3-24+2 
X3 = —12yı:18 ee SSS 
3 
Remove the equation (Fi,3). 
Consider f: = 2, t: = 2; go back to step 3. 


2x4-1 
Xı = — İX4t 6y13—-4 + 


= 
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Xat2yi3 
2322 yatai einen 
4 

X3 = —11yı3*8 (E2,3) 
2X4 + 3yzı - 1 = 0 (F2,1) 
Xa + 2yis — 4y22 = 0 (F2,2) 
Yes. 

Yes. İrə) = 1. We obtain x4 = -2yi3 + 4yz22, 
relation (T2). Substituting it in the others we 
get: 

—4Ayist8yoo-1l 
Xı = -28yas2 ot 20y13-4 + (E2,1) ' 
—3 

X2 = Y22 + yist5 (E2,2) ' 
X3 = = 11yı3t8 (E2,3) ' 


Remove the equation (F22). 





Consider f: = 3, t: = 
Xı = -22yis - 30y22-4 + 
X2 — yis + y22t5 
X3 = -llyış +8 


2yis + 2y22 + 3y31 — 1 = 


3 and go back to step 3. 


0 


—2yı3t2y22-1 





-3 


(Es,ı) 


(E3,2) 


(E3,3) 


(F3,1) 


10. 
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a = 1, and in = 1, jı = 31, j2 = 22 and r= 1. 
(A) yoz 21. —*Mö3ip relation (Hı). 


(B) Substituting it in the others we get: 





—2y13t221-2y31-1 





Xı = —-22yis -30Z: + 30y31 4 + (Es,1) 
=3 

X2 = yis + Zi — yait5 (E3,2) 

X3 = -llyış +8 (Es,3) 

2Y13 + 22, + Yap = 1 = 0 (F3,1) 

Xa = -2yis + Azı - 4yi3 (T2) 

(C) Consider h: = 2 


(B) ya = 1 — 2yi3 - 221, relation (Gi). 





Substituting it in the others we get: 





xi = -40yi3 - 922, + 27 (Es;1) ' 
X2 = 3yi3 + 3z:4 4 (E3,2) ' 
X3 = -llyis * 8 LESS" 
X4 = 6yi3 + 12z, - 4 (DM 
Yo. = 2yi3 + 32, - 1 (Hi) 


Remove the equation (F3,1). 





Consider g: = 2, f: = 4 and go back to step 3. 


Xy = —40yı3 = 92721 
X2 = 3yis + 3Z1 
X3 — -llyış 


No. The general 


system is: 
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+ 27 (E4,1) 
+ 4 (E4,2) 
+ 8 (E4,3) 


integer solution of the initial 


xi = -40k: - 92k, + 27, from (Eq.ı) 
X2 = 3ki + 3k. + 4, from (Ea,2) 
X3 = -11k +. 8; from (E4,3) 
Xa = 6ki + 12k, - 4, from (T2)'' 
X5 = 3k x 2, from (Tı) 
where kı, k2 0 Z. 

Algorithm 4 

Input 

A linear system (1) with not all ai; = O. 
Output 


We decide on the possibility of an integer solution 


of this system. 


integer solution. 


Method 


If it is possible, we obtain its general 
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(A) Divide every equation i by the maximal co- 
divisor of the coefficients of the 
unknowns. If the quotient is not an 
integer for at least one it, then the 
system does not admit integer solutions. 
Stop. 

(B) If there is an inequality in the system, 
then it does not admit integer solutions. 

(C) In case of repetition, retain only one 
equation of that kind. 


(D) Remove all the equations which are 


identities. 
Calculate a = min {|ai;|, ai; .. 0) and determine 
i,j 


the indices io, jo for which this minimum can be 


obtained. (If there are more variants, choose 





one, at random.) 

If a. 1, go on to step 6. 

If a = 1, then: 

(A) Calculate the value of the variable Xi, 
from the eguation iç; write this relation 
(Vv). 

(B) Substitute this relation (where possible) 


in all the equations of the system as well 


(C) 
(D) 
Does 


(A) 


(B) 


Write 


Write 


(Hy) . 
in al 


the r 
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as in the relations (Vy-1), (Hn) for all v 


and h. 





Remove the equation io from the system. 
Consider v: = vtl. 

at least one equation exist in the system? 
If it does not, write the general integer 
solution of the system substituting 

kı, ko, ... for all the variables from the 
right term of each expression representing 
the value of the initial unknowns of the 
system. 

If it does, considering the new data, go 


back to step 2. 


all ai j 5 Jo and bi, under the form: 


aij, Mit Yi, with [rs < EI ; 


da di + fior with İr: | « lass | 


Xj — x Ö qi j Xj + qi + T; relation 
dade © 2 
Substitute (where possible) this relation 


l the equations of the system as well as in 


elations (V,), (Hn) for all v and h. 
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8. Consider 
Xj = tu, h: = h + 1, 
0 
8i "EC Ti jr j "Sr 
deck ediş bDi Leş 
070 ES 070 0 0 


and go back to step 2. 


The Correctness of Algorithm 4 





This algorithm extends the algorithm from [4] 
(integer solutions of equations to integer solutions of 


linear systems). The algorithm was strictly demonstrated 





in our previous article; the present one introduces a new 


cycle--having as cycling variable the number of equations 





of the system--the rest remaining unchanged; hence, the 





correctness of algorithm 4 is obvious. 


Discussion 
Ls The counter variables h and v count the 
relations H and V, respectively, differentiating 
them (to enable the substitutions); 
2s Step 2 (A + B) + (C)) is trivial and is meant to 
simplify the calculus (as algorithm 2); 


3. Substep 5(A) has aesthetic function (as all the 
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algorithms described). Everything else has been 
proven in the previous chapters (see [4], [5], 


and algorithm 2). 


Example 4. Let us use algorithm 4 to calculate the 


integer solution of the following linear system: 


Solution 


3X1 = 7X3 + 6X4 = -2 


4X1 + 3X2 + 6X4 5X5 = 19 


a=3 and i = 1, j= 1 


3 m 1. Go on to step 6. 


So, 
= mS iy ck 
6 =3 ! 2 + 0 
-2 = 3 ! 0 — 2 
Xı = 3x3 — 2x4, + tı, relation (Hi). Substituting 


it in the second equation we get: 


Atı + 3x2; + 12xə — x4, — 5x5 = 19 


Xi: = ti, h: = 2, aio: = 0, ai3: = +2, aia: = Ü, 
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Go back to step 2. 


The equivalent system was written: 


+ 3ti + 3X3 = -2 


Atı + 3x2 + 12xə — x4, - 5x5 = 19 


(A) Then: x, = 4t, + 3x; + 12xə - 5x5 - 19, 
relation (Vi). 


(B) Substituting it in (Hi), we get: 


xi = - 7ti — 6x; - 21xs + 10xs + 38, (Hi) 


(C) Remove the second equation of the system. 





(D) Consider: v: = 2. 


Yes. Go back to step 2. 


The equation + 3t; + 2x3; = -2 is left. 
a=2andi=1, j= 3 

2 . 2, go to step 6. 

+3=4+2!2- 1 

— 2-12 (-1) + 0 

X3 = — 2t: + te — 1, relation (Hə). 


Substituting it in (Hı)", (Vi), we get: 
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Xi = 35tı - 6x2 — 21t> + 10x5 + 59 (Hi) "". 
Xa = — 20tı + 3X2 + 12t>2 = 5X5 d 31 (Vi) ". 
8. X3: = to, h: = 3, ayı: = -1, ais: = +2, bi: = 0 
(the others being all = 0). Go back to step 2. 
2. The equation -5: + 2t? = 0 was obtained. 
3. a = 1, and i = 1, j= 1 
4 1-1 


(A) Then, ti = 2t2, relation (Və). 


(B) After substitution, we get: 


Xi = 49tç, >> 6X» + 10x; + 59 (Hi) wt My 
Xa = —28t> + 3x2 - 5x5 — 31 (Vi) ''; 
X3 = — 3t» = 1 (H2)' ; 


(C) Remove the first equation from the system. 
(D v: = 3 
ss No. The general integer solution of the initial 


system is: 


Xı = 49ki — 6k; + 10k; + 59 
X2 = kə 

X3 = — 3ki = GL 
x4, = —28k: + 3k, - 5k3 - 31 
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where (ki, ks, k3) 0 Z’. 


Stop. 


Algorithm 5 





Input 


A linear system (1). 


Output 


We decide on the possibility of a integer solution of 
this system. If it is possible, we obtain its general 


integer solution. 


Method 


Die 


We solve the common system in R". If it does 
not have solutions in R", then it does not have 
solutions in Z” either. Stop. 

f=1, v=1, hzl 

Write the value of each main variable x; under 
the form: 


(Eri)i? Xi = Ó Qqij3X34 + qi + (Ó Ta TED) i 
j J 


with all Qiş, di, Tis, Ti, a. from Ly so that all 


İr::l < [8 [s |x: | < lal, “5. 0 (where all x;'s 


(Ye,i) 
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of the right term are integer variables: either 
from the secondary variables of the system or 


the new variables introduced with the 


algorithm). For all i, we write ri, / 5:. 
E 
(Ff, i)i: Ó Tiı:X: — Pij Y£,i + ri = 0, vhere 
b f 
J 
are auxiliary integer variables. Remove all the 


equations (Fr,ı) which are identities. 

Does it exist at least one equation (F¢,i)? If 
it does not, write the general integer solution 
of the system substituting ki, ko, ... for all 
the variables of the right member of each 


expression representing the value of the initial 
unknowns of the system. Stop. 


(A) Divide each equation (Fr,ı) by the maximal 
co-divisor of the coefficients of their 
unknowns. If the quotient is not an 
integer for at least one io, then the 
system does not admit integer solutions. 
Stop. 

(B) Simplify--as previously ((A)) all the 


fractions in the relations (Er,i)i. 


Calculate a = min (İrul, rij; .. 0), and determine 
3 


the indices io, jo for which this minimum is 
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obtained. 


If a . 1, go on to step 9. 


If a — 1, then: 
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(A) Calculate the value of the variable Xi, 
from the equation (Fg); write (V,) for 
this relation. 

(B) Substitute this relation (where possible) 
in the relations (Es), (Vwı), (Hn) for all 
i, v and h. 

(C) Remove the equation (Fri). 

(D) Consider v: = v + 1, f: = f + 1 and go back 
to step 3. 

9. Write all Ti, j .. Jo and ri under the form: 
fi $ EAST ! dig t fig; with Ix: al < àl; 

fi j ES Değe + ri with |e: | < Jal 

10. (A) Write See =- Ó 2552 + thy 
0 də 0 0 
relation (Hi). 

(B) Substitute this relation (where possible) 
in all the relations (Egi), (Fe,i), (Vv), 
(Hn-1) . 

(C) Consider h: = h + 1 and go back to step 4. 

The correctness of the algorithm is obvious. It 





consists of the first part of algorithm 3 and the end part 
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of algorithm 4. Then, steps 1-6 and their correctness 
were discussed in the case of algorithm 3. The situation 
is similar with steps 7-10. (After calculating the real 


solution in order to calculate the integer solution, we 





resorted to the procedure from 5 and algorithm 5 was 
obtained.) It means that all these insertions were proven 


previously. 


Example 5 
Using algorithm 5, let us obtain the general integer 


solution of the system: 


3X1 + 6X3 + 2X4 = 0 
4x2 = 2X3 = 7X5 = —1 
Solution 
1. Solving in R”, ve get: 
-—6X3-2X4 
Ri 
3 
2X%3+7Xs5-1 
KS E 
4 


10. 


(Er,ı) 


(E1,2) 


—2X4 
Xı = 2X3 "BS ———— 
3 
2X313X5—1 
X2 = X5 o ə ————— 
4 


—2 Xa kag 3yıı - 0 


2X3 + 3x5 — Ayı? — 1 = 0 


113 





2- wl 
322 ! 1+1 
-4 = 2 ! (-2) 
-l1=2 !0-1 
X3 = -X5 + 2yi2 + ti, relation (Hi). After 
substitution: 
—2X4 
(Ei,1)' Xi = 2X5 — 4yi2 — 2t1 + ———— 
3 
X5t4yiet2ti-1l 
(Ej,2)' X2 = Xs + 
4 
(Fi,2)' Xs + 2ti — 1 = 0 
Consider h: = 2 and go back to step 4. 
(Fii)' -2x, — 3yu = 0 


114 


(Era) r: 2ti + x — 1-20 


Yes 


a = 1 and io = 2, Je = 5 
(A) X5 = —2t: + 1, relation (Vi) 


(B) Substituting it, we get: 


—2Xa4 
(Eji)'': xi = -6ti + 2 — 4yı2 + ——— 
3 
(Ej2)'': x2 = -2ti + 1 + yo 
(H)! d$ x3 = 3ti— 1 + 2yi 


(C) Remove the equation (Fi,2). 





(D) Consider v = 2, f = 2 and go back to step 


3i. 
—2Xa4 
(Es; 1) : Xi = —6t. — 4yı? + 2 + ———— 
3 
(Ez,z) : x: = —2ti + yı + 1 
(F2,1) : -2x4 — 3y21 = 0 


a= 2 and ip = 1, jo = 4 
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(A) Xa = —yz: + te, relation (Hə). 


(B) After substitution, we get: 


2y»1-2to 
(E2,1) ' : Xı = —6tı — Ayı> + 2 +. —— SSS 
3 
(F2,1) ' : -ya - 2t2 = 0 
Consider h: = 3 and go back to step 4. 
(F2,1) ' : -ya - 2t2 = 0 
Yes 
a= 1 and io = 1, jo = 21 (two, one). 
(A) Vol = —2t>, relation (v2) i 
(B) After substitution, we get: 
(E21) '' : x1 = —6tı — Ayı, — 2t> + 2 


(H2) " : Xa = Gto 


(C) Remove the equation (F2,1). 





(D) Consider v = 3, f = 3 and go back to step 


(Es,1) : Xı = —6tı — diyi, — 2t2 + 2 


(Es,2) : X2 = -2ti + Yız + 1 


is: 


-6ki 
-2ki 


3k; 


-2ki 


where 
Stop. 


Note 1. 


(ki, 


+ 2k; 


k2, 


Algorithms 3, 
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The general integer solution of the system 


2k3 + 2 , from (Esa); 
+ 1 , from (E3,2); 
—13 5 from (Hı)", 
3k3, from (H2)'; 
+I y from (Vi); 
0 Z. 


4 and 5 can be applied in the 


calculation of the integer solution of a linear equation. 


Note 2. 


The algorithms, 


because of their form, are 


easy to introduce in the computer. 


Note 3. 


algorithm to use. 


It is up to the reader to decide on the 


Good luck! 
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